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CROSS REFERENCE TO RELATED APPLICATIONS 
5 This application is based upon and claims the benefit of priority from the Japanese 

Patent applications No. P2000-292770, filed on September 26 7 2000; the entire contents of 
which are incorporated herein by reference. 
BACKGROUND OF THE INVENTION 

1. Field of the Invention 

1U The present invention concerns the circuit design of a semiconductor device Such as 

LSI and is particularly for application in the estimation of a fluctuation of the circuit 
characteristic, as generated by process fluctuation, 

2. Description of the Related An 

s p5 Securing a high yield rate from the early stages of mass production is very important 

for the semiconductor devices, especially, LSI products. The reason is that the high yield 
rate leads to reduction of cost and turnaround time (TAT). This yield rate rises or falls due to 
fluctuation of the circuit characteristics. The factor which causes circuit characteristic 
fluctuation is manufacturing process fluctuation. If the manufacturing process did not 
lO fluctuate at all. the circuit characteristic also would not fluctuate, and if the circuit operation 
V could be thus assured, the yieJd rate of the LSI product could always kept at 100%. 
K However, in practice, it is impossible to remove completely manufacturing process 
fluctuation. 

In order to obtain a high yield rate, it is necessary to estimate the fluctuation of the 
25 circuit characteristic due to the fluctuation of the manufacturing process with high accuracy 
and optimize the circuit based on that estimation result- 
However, the estimated device and circuit characteristics sometimes do not coincide 
with the actual, measured values of the semiconductor device. That is to say, estimated 
values which emerge in simulation, may not always correspond to those found in practice. 
30 SUMMARY OF THE INVENTION 

The method for simulating the circuit characteristic of the semiconductor device 
according to the embodiment of the present invention is the obtaining a first characteristic 
value which most fluctuates a characteristic of an element composing the semiconductor 
device according to a fluctuation uf a parameter of the element; the determining a width of the 
35 fluctuation of the parameter which mutches u second characteristic value of the worst case of 
the characteristic of the element with the first characteristic value; and the determining a third 



characteristic value of the worst case of the circuit characteristic of the semiconductor device 

based on the width of the fluctuation. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a flow chart of a simulation method for the circuit characteristic of a 
Semiconductor device; 

Fig. 2 is a diagram for explaining a device structure parameter, which undergoes 
process fluctuation in the MOSFET (semiconductor device); 

Fig. 3 is a table for specifying the process fluctuation of the device structure 
parameter in each case of the worst-case analysis; 

Fig. 4 is the results of a calculation of the drive current (Idr) of CMOSFET (n type 
FET and p type FET) according to worst-case analysis in which the fluctuation width (K 
value) of the device structure parameter is 3; 

Fig. 5 is a structural diagram of a simulator for the circuit characteristic of the 
semiconductor device; 

Fig. 6 is a flow chart of the simulation method for the circuit characteristic of the 
semiconductor device; 

Fig, 7 is the results of a calculation of the drive current (Idr) of CMOSFET by 
according to worst-case analysis in which the fluctuation width (K value) of the device 
structure parameter is 1 ; 

Fig. 8 is the results of a calculation of the drive current (Idr) of CMOSFET according 
to worst-case analysis in which the fluctuation width (K value) of the device structure 
parameter is 2; 

Fig. 9 is a diagram showing the influence of the device structure parameter on the 
drive current (Idr) of the CMOSFET calculated by worst-case analysis; 

Fig. 10 is a diagram showing the drive current (idr) of the CMOSFET where the 
device structure parameter is fluctuated in a range from -3 o to +3 o ; 

Figs. 11 and 12 are diagrams for explaining optimization method for the fluctuation 
width (K value), which matches the results of the worst-case analysis with the results of 
Monte Carlo analysis; 

Fig. 13 is the results of a calculation of the drive current (Idr) of the CMOSFET 
according to the worst-case analysis in which the K value of the device structure parameter ts 
an optimized 2.375; 

Fig. 14 shows the circuit characteristic (access time) of each Case calculated based on 
the woist-case analysis of the device structure parameter where the fluctuation amount width 
(K value) is constant; 

Fig. 15 Is a diagram showing the circuit characteristics (access time) of ease* I and 



16, which are worst cases calculated by changing the fluctuation amount width (K value) of 
the device structure parameter; and 

Fig. 16 is a flow chan of the manufacturing method of the semiconductor device 
DETAILED DESCRIPTION OF THE EMBODIMENTS 

Various embodiments of the present invention will be described with reference to the 
accompanying drawings. It is to be noted that the same or similar reference numeral* arc 
applied to the same or similar parts and elements throughout the drawings, and the description 
of the same or similar parts and elements will be omitted or simplified. 

The present inventor sought for ihe cause of the estimated device characteristic and 
circuit characteristic not coinciding with actually measured values of the semiconductor 
device. That is, the reason for estimated device characteristic and circuit characteristic that 
are possible in simulation not being found in practice will be clarified. 

The present inventor considers that the reason these values do not coincide lies in the 
accuracy of worst-case analysis on the device characteristic being low, and compares these 
results with Monte Carlo analysis, which is capable of calculating distribution characteristics 
due to process fluctuation of the device characteristic at high accuracy. The "worst-ca-.e" 
mentioned here refers to a case in which maximum fluctuation in manufacturing process 
occurs. 

Fig. 1 shows a flow chart of the simulation method for the circuit characteristic of a 
semiconductor device used to search for the cause. First, us shown m step Sll, a fluctuation 
width of manufacturing process of the semiconductor device is input. Next, in step Si 2, 
worst-case analysis on the process flucTuaiion is carried out. In this analysis, the fluctuation 
width of the electric characteristic of each device composing the semiconductor device is 
obtained. Next, in step S13 7 the worst-case analysis on an electric circuit in the 
semiconductor is carried out based on the fluctuation width of the electric characteristic of the 
device. From this analysis, the circuit characteristic of the semiconductor device is obtained. 

Because main stream semiconductor devices today, such as LSI 7 are complementary 
metal oxide semiconductor field effect transistors (CMOSFET), these CMOSFET will be 
described here. The CMOSFET is a technology for constructing a circuit by 
complementarily combining an n-type MOS field effect transistor (nFET) with a p-typc MOS 
field effect transistor (pFET). Therefore, the LSI circuil is constructed from the n-type FET 
and the p-iypc FET as its basic devices, so that they are integrated with each other. As 
evident from thi* fact, the cause of the circuit characteristic fluctuation is the device 
characteristic (electric characteristic of the n type FET and the p type FET) fluctuation, 
Further, the cause of the device characteristic fluctuation is that process fluutuatiun* 
(fluctuation of process parameter) are produced in the p type FET and the p type FET during 



manufacturing. The device characteristic mentioned here refers to the electric characteristics 
of the n type FET and the p type FET, for example, the threshold voltage or saturated current 
value. The circuit characteristics refer to so-called the circuit specification, for example, 
propagation delay time of circuit, access time, frequency, power consumption and the like 

The n type FET and p type FET of Fig. 2 each contain a silicon substrate 31, The 
silicon substrate 31 includes a source region 33 and a drain region 34, which is located 
separate from the source region 33 but so that both of them include the surface of the silicon 
substrate 31, Further, the silicon substrate 31 includes a channel region 32 between the 
regions 33 and 34. A gate insulation film 35 is provided on the silicon substrate 31. A gate 
electrode 36 is provided above the channel region 32 on the gate insulation film 35, 

There are three important device structure parameters for determining the device 
characteristic in the n type FET and the p type FET shown in Fig. 2, namely, gate length (L), 
gate oxide film thickness (Tox) and channel impurity concentration (Nch), Then, in step Sll, 
us the fluctuation width of process, fluctuation width (AL) of gate length (L), fluctuation 
width (A Tox) of gate oxide film thickness (Tox) and fluctuation width (A Nch) of channel 
impurity concentration (Nch) are input. 

Next, in step S12, worst-case analysis of the process fluctuation is carried out. 
Through the worst-case analysis, circuit characteristic fluctuation caused by the process 
fluctuation can be estimated. The basic concept of the worst-case analysis is that the worst 
case of the process analysis provides the worst case ot the device characteristics. And the 
worst case of the device characteristics provides the worst case of the circuit characteristic. 
Fluctuations in the device characteristics and circuit characteristics are estimated. First, the 
worse case of the process fluctuation is obtained, 

The worst case of the process fluctuation can be generated by combining the 
minimum value and the maximum value of each parameter P(L ? Tox, Nch), Generally, the 
fluctuation of the process parameter indicates Gaussian distribution horizontally symmetrical 
to a nominal value (a process parameter value when the fluctuation amount is zero). If the 
nominal value in the Gaussian distribution concerning process parameter P is PU, standard 
deviation is o , according to the principle of 3 0 in statistics, and the probability that the 
fluctuation of the process parameter P falls within P0±3 o is 99.7%, Therefore, it may be 
considered that the minimum value of the process parameter is PU-3 o and ihc maximum 
value is PO+3 O . Then, in the worst case analysis, fluctuation amounts +3 o and -3 o are 
used as the maximum value and minimum value of each process fluctuation. That is, if the 
nominal value and standard deviation in Gaussian distribution of each of the gate length (L) 5 
gate oxide film thickness (Tox) and channel impurity concentration (Nch) are (LU, a L ), (TO, 
(J r )> (NU, a N ) 3 the minimum value and maximum value of each process parameter are (Lt) - 



3 o L , LO + 3 o L ), (TO - 3 o T , TO + 3 o T ), (NO - 3 o N , NO + 3 o 

The two process fluctuation amounts A L and A Tox are the same as the fluctuation 
amount in the n type FET and p-type FET in the same chip. However, the fluciuaiion 
amount ANch has a completely independent fluctuation amount to each of the other n-type 
5 FET and the p-type FET despite being part of the same chip. That is 7 in the CMOS circuit, 
channel impurity concentration fluctuation of the n type FET (ANch (n type FET)) and 
channel impurity concentration fluctuation of the p type FET ( A Nch (p type FET)) need to be 
handled independently of each other. Thus, the process fluctuations which should be 
considered in order to estimate the circuit characteristic fluctuation comprise four types of the 
10 fluctuations, namely, AL, A Tox, A Nch (nFET) and A Nch (pFET). 

In the worst-case analysis, as shown in Fig. 3 ? the device characteristic and the circuit 
characteristic are calculated for a total of 16 cases (2 4 ~ 16) in which the maximum value (+) 
Q and minimum value (-) of each process fluctuation of four types. Then, the worst case is 
« selected from the 16 cases. More specifically, a process fluctuation shown in Fig. 3 (case: 1 
\§i — 16) is given to nominal spice model (input parameter set for use in circuii simulation) 
j~ corresponding to nominal condition. That is, by adding a fluctuation amount AP to the 
% nominal spice parameter P0, namely P = P0 + A P 7 a new corner spice model P is generated, 
ry Consequently, 16 comer spice models are generated, and then simulations of the device 

characteristics and circuit characteristics can be executed for each of the 16. 
2tt A result of the worst-case analysis on the device characteristic is compared with 

W Monte Carlo analysis capable of calculating the distribution characteristic of the device 
p characteristic by process deviation at a high accuracy. 

IM* Fig. 4 indicates wiih a square where each case shown in Fig. 3 is located on an Idr 

distribution diagram of the device characteristic. Its abscissa (x-) axis indicates the drive 

25 current (Idr (n type FET)) of the n type FET while the ordinate (y-) axis indicates the absolute 
value (|Idr (p type FET)|) of the drive current of the pFET. Fig. 4 shows a case where all + 
in Fig. 3 are fluctuation width of +3 o while all - indicates fluctuation width of -3 O . 

As a reference for the worst-case analysis, the strictly defined Idr distribution 
diagram obtained from a circuit simulation by Monte Carlo analysis is indicated by the ellipse 

30 formed by connecting symbols* X. According to Monte Carlo analysis, the Idr value can not 
obiain any other value than values within this ellipse. Tln= ellipse of Xs rising in a righiward 
direction as a periphery to the distribution diagram is obtained by the four process fluctuations 
(AL, A Tox, ANch (nfet)) each being provided with independent Gaussian distribution and 
then the Monte Carlo calculation being executed 2,OUU times by circuit simulation. Because 

35 the Idr distribution diagram indicates the ellipse rising in the rightward direction, it found 
that there is a strong correlation between Idr (nfut) and |ldr(pfct)|. Further, crubs-Iiice arrows 



14, 17 located slightly lower and left of the center of the ellipse in each diagram denote the 
nominal value (vertical arrow 17) of Idr (nlet) and the nominal value (lateral arrow 14) of |Idr 
(pfetj|- 

The reason why the worst-case analysis is executed although the fluctuation width of 
the circuit characteristic can be found just by executing Monte Carlo analysis is that the 
analysis time of Monte Carlo analysis is very long. Because a relatively small circuit 
(smaller than several thousand devices) dose not require so long a simulation time, it is 
permissible to obtain its fluctuation values directly by using Monte Carlo analysis on the 
circuit characteristic. However, in the case of a large-scale circuit (several thousand - 
several ten thousand devices), because a single simulation time may sometimes take several 
hours, application of the Monte Carlo analysis tor example 2000 times is not practical in view 
of calculation efficiency. This amount of analysis lime would be fatal to the development of 
a semiconductor in which development speed is demanded, and therefore, Monte Carlo 
analysis cannot be applied. 

On comparison of both the analysis results, it is made evident that some of the 
squares which indicate the results of the wor>st-case analysis are largely deviated from the 
strictly defend elliptical Idr distribution diagram based on Monte Carlo analysis. The 
existence of these deflected squares is a fatal problem for the worst-case analysis. Thai is, 
by right the process fluctuation should have a fluctuation width of ± 3 o with resped to 
each process parameter. However, the probability ihat the absolute values of respective 
process parameters become 3 O at the same time is very low, and in practice such a case 
hardly ever occurs. Thus, in the worst-case analysis shown in Fig. 4, an excessive 
(over-estimated) estimation result is obtained. This is because worst-case analysis does not 
fit practical situations. 

Finally, in step Sl3 3 from the excessive estimation result, a corresponding comer 
spice model is obtained, and then using thai corner spice model the worst case of the circuit 
characteristic is estimated. In this case, the excessive estimation result of the device 
characteristic produces an excessive estimation result of the circuit characteristic. The 
estimation accuracy of the circuit characteristic is not assured at all. This is a problem with 
the worst-case analysis. That is, the worst-cuse analysis is not capable of estimating al high 
accuracy the device characteristic due to a process fluctuation and u fluctuation of the circuit 
characteristic. 

(Organization ot a simulator for the circuit characteristics of a semiconductor device) 

Fig. 5 is an organization of a simulator tor the circuit characteristics uf a 
semiconductor device according to this embodiment. The Mmulutor for the circuit 



characteristic of the semiconductor encompasses a computing unit 1, recording units 6 to 9, 
and input/output unit 10, The computing unit 1 7 the recording units 6 to 9 and the 
input/output unit 10 are connected to each other through a bus 11. The computing unit 1 
includes a process fluctuation Monte Carlo analyzing unit 2, a worst-case analyzing unit 3 for 
5 process fluctuation to arbitrary K-value 7 a computing unit 4 for K-value corresponding to 
maximum fluctuation width, a worst-case analyzing unit 5 for circuit characteristic, a circuit 
characteristic determining unit 41, a process parameter selecting unit 42, a specification 
setting unit 43 tor device characteristic and circuit characteristic and a resetting unit 44 for 
nominal value and fluctuation width of device structure parameter- The recording units 6 to 
10 9 include a circuit information recording unit 6, a process fluctuation width recording unit 7, a 
device characteristic fluctuation width recording unit S 5 and a circuit characteristic fluctuation 
width recording unit 9. 

yj (Simulation method for the circuit characteristic of semiconductor device) 

If Fig. 6 is a flow chart of a simulation method using the simulator about the circuit 

1g characteristic of the semiconductor device according to this embodiment. 

S First, in step Si, wiring information of the semiconductor device and circuit 

information about circuit size such as iho quantity of the FETs are input through the 
b input/output unit 10 and recorded in the circuit information recording unit 6 as circuit 
W information. Further, as proven frequency distribution values, the parameter of the device 
m structure such as the FET is input through the input/output unit 10. If this frequency 
iP distribution is recognized to be Gaussian distribution, the nominal value and the standard 
£2 deviation G are input. This frequency distribution is recorded in the process fluctuation 
width recording unit 7. For the process fluctuation width, fluctuation width (AL) of gate 
25 length (L), fluctuation width (ATox) of gate oxide film thickness, and fluctuation width (A 
Nch) of channel impurity concentration (Nch) are used. 

Meanwhile, the frequency distribution of the device structure parameter about the 
gate length (L)> gate oxide film thickness (Tox) and channel impurity concentration (Nch) are 
regarded us Gaussian distribution- Thus, the fluctuation width from each nominal value of 
3D the device structure parameter can be expressed such that the standard deviation o is a unit. 
The fluctuation width expressed such that the standard deviation o is a unit is called K-value, 
Then, if the fluctuation width of each device structure parameter is converted to a K-value, the 
procedure after step S5 can be easily executed. Thus, each fluctuation width of plural device 
structure parameters can be expressed as a K-value, Meanwhile, the K-value may be 
35 different depending on the device structure parameter. For example, the K-value may be ^>et 
to be the same value or a different value where it increases from tiie nominal value or 

7 



decreases from the nominal value, corresponding to the characteristic of the process, 
likewise, the same K-value may be set up as the gate length (L), gate oxide film thickness 
(Tox) and channel impurity concentration (Nch) or different K-values may be set up. If 
different K-values are set up, they are handled so that they can be distinguished from each 
other. 

Next, in step S2, it is determined whether or not Monte Carlo analysis of the circuit 
characteristic of the semiconductor device is appropriate. If possible, the Monte Carlo 
analysis on the circuit characteristic in step S3 is carried out. The is a criterion for 
determination is that time required for Monte Curio analysis on the circuit characteristic 
should not be too long. More specifically, the circuit size of the FET or the like composing 
ihe circuit acts as the criterion for determination. In step 3 3 the Monte Carlo analysis of" the 
circuit characteristics, and in step S4 7 the Monte Carlo analysis of process fluctuation on 
individual devices such as the FETs composing the circuit are carried out. In step S4, the 
Monte Carlo analysis of the process fluctuation is carried out in the Monte Carlo analyzing 
unit 2 for process fluctuation shown in Fig. 5. In this analysis, the fluctuation of the circuit 
characteristic is obtained, and then the simulation method is complete. 

If it is determined in step S2 that the Monte Carlo analysis of the circuit 
characteristic is inappropriate, steps S5, S6 are carried out. Although the time necessary for 
this simulation can be reduced by executing both steps in parallel, it is possible to execute 
either of them ahead of the other. 

In step S5, the worst-case analysis of the process fluctuation with K-value as a 
variable is carried out on devices by the worst case analyzing unit 3 tor process fluctuation 
corresponding to arbitrary K value ot Fig. 5. The circuit of the semiconductor uses plural 
devices such as the FETs. Because only the same result can be obtained even if process 
fluctuation analysis is carried out for the same structure device, the analysis is carried out 
once for each different structure. For example, if the nFETs and nFETs in the circuit are all 
of the same structure, the analysis on the nFET and pFET only has to be carried out once for 
each. 

Figs, 7 and 8 show the results of a calculation of driving current Idr corresponding to 
process fluctuation by changing the fluctuation amount width (changing the K-valuu) by 
carrying out worst-case analysis of the nFET and pFET Fig. 7 shows a case where the 
K-value is ±1, so that the fluctuation width of gate length (AL), the fluctuation width of 
gate oxide thickness (A Tox), and channel impurity concentration fluctuation width (A Nch) 
are set to the same K value ±1. Further, the fluctuation width from the nominal value in a 
positive direction and in a negative direction is set as equal for the K value +1 and the K value 
- 1, On the uther hand, Fig. 8 indicates a case where the K value is ±2. Fig. 4 described 

8 



previously indicates a case where the K value is ± 3. 

In Fig. 3, the location of the respective cases 1 to 16 on the Idr distribution diagram 
as device characteristics are indicated by the corresponding squares. The case Nos, 1 to 16 
correspond to Cl to C16 in Figs, 4, 7 and 8. The abscissa (x) axis in each of Figs. 4, 7 and 8 
5 indicates an absolute value (Idr(nFET» of a drive current corresponding to nFET, while the 
ordinate (y) axis thereof indicates an absolute value (|Idr(pFET)|) of the drive current 
corresponding to pFET. Because the intersection of arrows 14 and 17 provides a drive 
current with no process fluctuation, the square most distanced from this intersection provides 
the worst case. From this 7 it is found that when the K value is ±1, case Cl and C16 are the 
10 worst cases. Even when the K value is ±2 in Fig. 8 and ±3 in Fig. 4, cases Cl and C16 
are the worst cases. 

The fluctuation width W of the device characteristic refers to the distance between 
Cl thc inlerseCti on of arrows 14 and 17 and case 1 or 16, or the distance on the graph plane of 
; C Figs. 4, 7 and 8, The driving current of the worst cases Cl, Cl 6 and the fluctuation width W 

]Jg of rite device characteristic are recorded in the device characteristic fluctuation width 
JS recording unit 8 of Fig. 5. These are compared to Figs. 4, 7 and 8. In case Cl, it is found 
^ thai Idr(nFET) and |Idr(pFET)| decrease as the K value changes to -1, -2, -3. About the ca*e 
m C16, it b abo found that Idr(nFET) and |ldr(pFET)| increase as the K value changes to 1, 2, 3. 
^ Fig. 9 is a diagram showing the relationship between the driving current of each case, 

M as calculated by worst case analysis of the process fluctuation in CMOSFET and the 

"Saw' 

ifM fluctuation amount of each process fluctuation. The disposition position of each case is a 
point which can be controlled depending on the size of a process fluctuation (AL, A Tox 3 A 
y ; Nch (nFET), A Nch(pFET)). The 16 cases Cl to C16 arc divided into four groups Gl to G4. 
For example, the group G4 is composed of the four ca.se* C13 to C16 and the disposition 

25 positions of these four cases C13 to C16 allow a rectangle to be formed with these disposition 
positions as corners. The positions of groups Gl to G4 are determined by the sum (AL+A 
Tox) of the fluctuation amounts of AL and A Tox. The degree of spreading of the four 
cases in each of the four groups in the direction of Idr(nFET) is determined by twice A Nch 
(nFET). The degree of spreading in the direction ot |Idr(pFET)[ is determined by twice A 

30 Nch(pFET). From these, it is found that the fluctuation amount of Idr by the process 
fluciuation of AL and A Tox is larger than the fluctuation amount of Idr by the process 
fluctuation of A Nch (nFET) and A Nch (pFET). The reason enabling The comparison of 
the size of the fluctuation amount of the device characteristic in relation 10 the fluctuation 
amount is the transformation of the fluctuation amount is to a. K-value, which is a standard 

35 vaJue. 

As for the CMOSFET, it can be easily found that the worst cases are Cl, C16 from 

9 



such a graph as Fig, 7, If the worst case is automatically determined, it is necessary to 
calculate the distance from an Intersection between the nominal values 14 and 17 for all the 
cases Cl to C16, so as to select the maximum value for this distance. However, the 
distances for ail the cases do not always have to be calculated. That is 7 in addition to the 
distance from the intersection between the nominal value of the case Cl and that of the case 
Cl 6 7 1/2 a distance from the position of the case C6 up to the position of Cll, 1/2 a difference 
of the ldr(nF6T) values between the case C6 and case Cll are to be chosen for the maximum 
value. The same thing can be said for case C7 7 CIO. Thus, it is possible to reduce the 
number of choices. 

Fig. 10 is a diagram showing the relationship between each process fluctuation and 
fluctuation amount of the driving current in the CMOSFEX This shows the values of Idr 
which change when each process fluctuation (AL 7 ATox, A Nch(nFET) 7 A Nch(pF£T» is 
fluctuated to —3, -% -1, 0 ? -t-l 5 +2 7 +3 by the K value. For example, when fluctuating the A 
L, other process fluctuation than A L is set to zero. The larger the absolute value of the K 
value, ihe farther the value of Idr upon the process fluctuation moves from the nominal value. 
It is found that the fluctuation of AL contributes to the fluctuation of Idr most largely, 
followed by the ATox, The A Nch(nFET) and AJNch(pFET) contribute to fluctuation 
only of Idr(nFET) and only Idr(pFET). Although the fluctuations of AL and ATox 
contribute to the fluctuations of Idr(nFET) and Idr(pFET), they hardly contribute to the 
fluctuations of A Nch(nFET) and A Nch(pFET). Consequently, by obtaining the 
fluctuations of ldr(nFET) and ldr(pFET) without considering . the fluctuations of A 
Nch(nF£T) and A Nch(pFET), the fluctuations of the device characteristic and circuit 
characteristic can be estimated more quickly while maintaining high accuracy. 

In step S6 of Fig. 6, Monte Carlo analysis of process fluctuation is carried out on a 
device having the same structure us the device subjected to analysis in step S5, by the Monte 
Carlo analyzing unit 2 for the process fluctuation. In Figs. 4, 7 and 8, the same results 
Calculated by Monte Carlo analysis on a driving current corresponding to the process 
fluctuation in the CMOSFET are expressed with a cross symbol. The reason the 
representation is repeated in the respective diagrams is as a reference for the worst-case 
analysis. 

The strict Idr distribution diagram obtained from Monte Carlo circuit simulation is 
indicated by an ellipse drawn with cross symbols. That is, 54 cross symbols disposed in the 
shape ot an ellipse rising in the nghtward direction of each diagram are located on the 
outermost periphery of a distribution diagram oi 2,()(X) cross symbols obtained by providing 
four kinds ot process fluctuations (AL, A To% 7 A Nch(nfet), A Nch(pfet)) each wiih 
independent Gaussian distributions, and then executing a Monte Carlo calculation 2,OOU times 
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based on a circuit simulation. 

Because the Idr distribution diagram indicates an ellipse rising in the rightward 

direction, it is evident that there is a strong correlation between Idr (nFET) and |Idr (pFET)|. 

The arrows 14 and 17 located slightly below and leftward of the center of the ellipse of each 
5 diagram denote a nominal value (vertical arrow 17) of Idr (nFET) and a nominal value (lateral 

arrow 14) ot (Idr (pFET)J. Meanwhile, the arrows 16 and 18 indicate a minimum value and a 

maximum value of |Idr (pFET)|. 

Because the intersection of the arrows 14 and 17 provides a drive current where there 

is no process fluctuation, the cross symbol most distanced from this intersection provides the 
10 maximum fluctuation width M of the device characteristic. Consequently, it is evident that 

the cross symbols, which provide the maximum fluctuation width M are cross symbols SS, FF. 

The maximum fluctuation width M of the device characterisxic refers to the distance on the 
Q graph plane between the cross symbols SS and FF. 

In Figs. 4, 7 and 8, the disposition positions of the cases Cl 5 C2 which are the worst 
S cases with respect to the disposition positions of the cross symbols SS, FF change depending 
^ on the K value. If the distance from the intersection of the arrows 14 and 17 is regarded a> a 
"5 reference point, when the K value is ±1 or ±2, the distance from cases Cl and C16 to the 
t& intersection is shorter than the distance from the cross symbols SS, FF to the intersection. 
% On the other hand, if the K value is ±3, the distance from cases Cl, C16 to the intersection 
20 is longer than the distance from cross symbols SS, FF to the intersection. Consequently, by 
r t setting the K value at a specific value located between -t-2 and +3 and a specific value iucated 

between -2 and - 3, it is possible to match or dispose the disposition positions of the cases Cl 3 

C2 7 which are the worst cases, with the disposition positions of ihc cross symbols SS, FF, or 

in the vicinity thereof. 

25 Then, the K value is optimized »o that the position of the worst case coincides with 

the edge of a proper Idr distribution diagram (elliptic). The optimization ot the K value 
mentioned here denotes optimizing the value of a parameter K when the maximum fluctuation 
amount is expressed as (maximum fluctuation amount) = ±K tJ. It does not denote 
optimizing a itself. As shown in step S7 of Fig, 6, the fluctuation width W in the worst case 

30 is brought near the maximum fluctuation width M in Monte Carlo analysis by changing the K 
value. Thai is, a K value which minimizes the value of a function F is obtained. For this, 
the fact that positions on the graph of cases Cl to C16 in the worst case analysis can be 
controlled by changing the size of the process fluctuations (AL ; ATox, A Nch(nFET), A 
Nch (pFET))(by changing the K value) is employed. 

35 More specifically, the positions of the cases Cl, C16 ? which are the worst vases, on 

the Idr distribution diagram are matched wirh the positions of the cross symbol SS (xOiowcr, 



yOlower) and the cross symbol FF (xOupper, yOupper) oq the strictly defined Idr distribution 
diagram obtained by Monte Carlo analysis. Meanwhile, although the values of xOlower and 
xOupper arc values of ldr(nFET), the Idr(nFET) value may be standardized on a scale in 
which the distribution is in a unit with the nominal value as zero. The values of yOlower and 
yOupper may be standardized also. Thus, the above-described alignment can be carried out 
in circuit characteristics having different unit systems and the like. Then, function F defined 
in the Eq T (1) is introduced to carry out the matching, 

F=[(x-xO) 2 + (y-yO) 2 ) ,/2 - (1) 

In the Eq. x and y denote positions (Idr value), (x ? y) on the Idr distribution 
diagram of the case Cl or Cl6 7 which is the worst case. X and y are expressed as function 
of the K value. On the other hand, xO is general expression of xOlower and xOupper and yU 
is general expression of yOlower and yOupper. As for the function F, the function Fl 
expressed by the Eq. (2) denotes the distance between the case Cl and the cross symbol SS on 
the Idr distribution diagram. Likewise, the function Fl 6 expressed by the £q. (3) denotes the 
distance between the case C16 and the cross symbol FF on the Idr distribution diagram. 
Consequently, the parameter K value which minimizes the values of the functions Fl, F16 to 
carry out the matching is obtained. For simplification here, the K values of the fluctuation 
values of the device structure parameters of four kinds are moved in uniform. 

Fl = [(xl - xOlower) 2 + (yl - yOlower) 2 ] 172 — (2) 
F16 = L(xl6 - xOupper) 2 + (y 16 - y0upper) 2 ] 1/2 - * (3) 

With the abscissa axis as an absolute value of the K value and the ordinate axib as the 
function F, the sum of the functions Fl and F16 and the sum of the functions Fl and F16 {Fl 
+ F16) are shown. Consequently, it is evident that the function Fl is minimized when the 
absolute value of the K value is 2.19. The function F16 is minimized when the absolute 
value of the K value is 2.43, The sum of the functions Fl and F16 (Fl + F16) is minimized 
when the absolute value of the K value is 2.375. The purpose for evaluation of the sum of 
the functions Fl and FIG (Fl + F16) is to find Out a single K value which reduces both the 
functions Fl and F16 of the Cl and C16, which are two poieb of the worst cases. Hereinafter, 
this K value of 2.375 is called value after optimization". As the K value after 
optimization, ii Is possible to use the average of the absolute value of the K. value which is 
minimized under the function Fl and the K value which is minimized under the function F16. 
Further, because the absolute value of the K value which is minimized under the function Fl 



is larger than the absolute value of the K value which is minimized under the function F16, it 
is possible to consider that case C16 is u worse case than case Cl, and employ the absolute 
value of the K value which is minimized under the function F16 as the absolute value 2.43 of 
the K value after optimization, 
5 Fig. 13 shows the results of obtaining the positions of the cases Cl to C16 on the ldr 

distribution diagram by worst-case analysis using a K value of 2.375 as the K value after 
optimization. In Fig, 13, the K values of the four kinds of process fluctuations are uniformly 
±2.375, as in Fig, 12. Using Monte Carlo circuit simulation, the edge of the strictly defined 
elliptic ldr distribution diagram can be well matched with the positions of the cases m the 
10 worst-case analysis. 

In summary, in step S6 of Fig. 6 7 the strict ldr distribution diagram ot device 
characteristic fluctuation is obtained by Monte Carlo circuit simulation (Monte Carlo device 
p simulation may also be used). The Monte Carlo analysis mentioned here is applied to the 
device characteristic of a single MOS transistor. Thus, particularly, in circuit simulation, a 
fall in computation efficiency, which is ordinarily problematic to Monte Carlo calculation, 
4£ does not have to be a cause of concern. Next, in step S5, the worst-case analysis is carried 
15 out on the device characteristic of the single MOS transistor. The analysis result is recorded 
fg in the fluctuation width recording unit 8 of the device characteristic ot Fig. 5, Then, by 
;|L matching the results of the worst case analysis with the resuhs of ihe Monte Carlo analysis as 
2§j shown in Fig. 12 using the K-value computing unit 4 corresponding to a maximum fluctuation 
iPW width of Fig. 5, optimization of the K value in the worst-case analysis is carried out, so that 
*m *he K value after optimization is obtained. 

In step S8 of Fig. 6, the processing step another advances by the worst-case analysis 
of the circuit characteristics of Fig. 5 being carried out by the worst-case analyzing unit 5, so 
25 as to obtain a fluctuation 25 in the circuit characteristics. This analysis result is recorded in 
the fluctuation width recording unit 9. In this case, using the " O value after optimization" 
already obtained from the fluctuation analysis of the device characteristic, the fluctuation 25 
of the circuit characteristic of a large scale circuit can be estimated quickly at high accuracy. 

The estimation result ot ihe fluctuation 25 of the circuit characteristic will be 
30 described with reference to Figs. 14 and 15. Fig. 14 shows the behavior of access time (t AC ) 
as the circuit characteristic. This figure shows value* of acccs* times to the cases Cl - C16 
when the K value m the worst case analysis is ± 1 7 ±2, ±:3. A maximum value (a 
horizontal arrow 21 of about 5fns]) of access time obtained by strict Monte Carlo analysis, a 
nominal value (a horizontal arrow 2<J of about 3,9[ns]) and a minimum value (a horizontal 
35 arrow 19 of about 3.4[ns]) are indicated on the ordinate axis. Consequently, it is evident that 
when the K value is ±3a, The maximum value of access time in ihe case Cl to C16 i.s 
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over-estimated, while the minimum value of access time in the case C16 is under-estimated. 
When the K value is ±1, ±2, the maximum fluctuation width from the nominal value is 
smaller than ihe minimum fluctuation width obtained by Monte Carlo analysis. This 
denotes that in order to estimate the fluctuation (maximum value and minimum value) of the 
5 access time as a circuit characteristic accurately according 10 the worst-case analysis, a K 
value or the like, whose absolute value is in a range of 2 or more to less than 3, must be used. 
This coincides with the results obtained in step S7. 

Fig. 15 indicates the behavior of the access time t A c (ordinate axis) when the absolute 
value (abscissa axis) of the K value is changed as a parameter for the cases CI and C16. 
10 From this 7 it is found that in order to estimate the maximum value of the strict access time in 
the case Cl, the absolute K value is set at 2.35. Likewise, it is found that in order to estimate 
the minimum value of the strict access time in the ca^e C16, the absolute K value is set at 2.4. 
q The absolute K value coincides with the absolute K value introduced according to the 
® behavior of the device characteristic in Fig. 1 2. That is, according to the procedure shown in 
fjpj| Fig. 6, the K value is optimized to a device characteristic fluctuation shown in Fig. 12 in steps 
«p S5 to S7, Next, in step S8, the worst-case analysis is carried out on circuit characteristic 
^ fluctuation using the "K value after optimization". Consequently, fluctuation 25 of the 
gj circuit characteristic originating from process fluctuation can be estimated quickly and at high 
^ accuracy, which cannot be achieved in worst-case analysis, 

|j§ By optimizing the circuit (optimization of the nominal value of the process parameter 

HI or optimization of the quantity of transistors) based on the result of fluctuation 25 of the 
JL" circuit characteristic obtained according to the simulation method for the circuit characteristic 
y* of a semiconductor device, the process window can be widened, thereby achieving a high 

yield rate in mass production of the LSI. 
25 As described above, it is possible to provide a simulator for the circuit characteristics 

of a semiconductor device and a simulation method capable of estimating fluctuation in the 

manufacturing process of the semiconductor device and fluctuation in the circuit 

characteristic quickly and at high accuracy, 

30 (Manufacturing method of semiconductor device) 

Next, a manufacturing method of the semiconductor device according to the 
simulation method for the circuit characteristic of the semiconductor device will be described. 
The manufacturing method of the semiconductor device occurs according to the flow of the 
simulation method for the circuit characteristic of the semiconductor device shown in Fig. 6, 

35 and then that of Fig. 16. In the description below, the procedure after the step S2I is 
executed using the obtained fluctuation 25 of the circuit characteristic, will be described, 
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First, in step S21 using circuit characteristic determining unit 41 of Fig, 5, it is 
determined whether the specification 26 of the circuit has been satisfied, even if the circuit 
characteristic of the semiconductor device fluctuates, according to the fluctuation 25 of the 
circuit characteristic. If the specification 26 is satisfied, the semiconductor device is 
manufactured without changing an initial phase nominal design value of the device structure 
parameter in step S25. Then, manufacturing of the semiconductor device is complete. 

When the device characteristic and circuit characteristic deflect from their 
Specifications 26, the following two methods are available. 
(Method 1) 

A possible value of the device structure parameter is an initial phase nominal value 
plus a fluctuation width originated from a process condition, For example, the gate length L 
is expressed in the expression ot L = LO + A L (LO; initial phase nominal design value, AL: 
fluctuation width). This fluctuation width A L originates from device manufacturing 
environment and accuracy/performance of manufacturing machine. Therefore, the 
fluctuation width A L is reduced by adjusting the manufacturing environment and the 
manufacturing machine. By this operation, disparities of the device characteristic and circuit 
characteristic are minimized, so that these characteristics can be maintained within the 
specification. However, reduction of the control width of the device structure parameter by 
adjustment of the manufacturing environment and manufacturing machine may be sometimes 
proved difficult. 
(Method 2) 

The initial phase nominal design value is changed. For example, in case of gate 
length L, the initial phase nominal design value LO is changed. Then, the device 
characteristic and the circuit characteristic are realigned so that they satisfy the specification. 
This method is an optimization method called design centering. 

Hereinafter, mainly the method 2 will be described. 

If the specification 26 is not satisfied, in step S22, the device structure parameter 
which in a cause for not satisfying the specification is selected by the device structure 
parameter selecting unit 42 of Fig. 5. For example, assume that as shown in Fig. 15, the 
upper limit value of the specification of access time tAC is 4.8 nanoseconds. Because the 
upper limit value of the access time tAC is Known to be 5.0 nanosecond from the fluctuation 
25 in the circuit characteristic, in step S21, it is determined that the specification is not 
satisfied. Because the specification Is not satisfied, the processing proceeds to step S22. In 
step S22, a device structure parameter in which the device characteristic and circuit 
characteristic change largely with respect to the K value as shown in Fig. 10 is selected, in 
Fig, 11), for example, the fluctuation of the gate length ( A L) is selected or fluctuations of A 



L and gate oxide film thickness ( A Tox) are selected. 

In step S23, the settings of the initial phase nominal value of the selected device 
structure parameter are realigned, so that the device characteristic and the circuit characteristic 
satisfy the specification, using the resetring unit 44 for the nominal value of the device 
structure parameter and the fluctuation width in Fig. 5. For example, the selected gate lengih 
L and the initial phase nominal value of the gate oxide film thickness Tox are reset to new 
values. 

Finally, in step S24, a semiconductor device is manufactured according 10 a newly set 
nominal design value of the device structure parameter. 

As described above, according to the present invention, it is possible to provide a 
manufacturing method of the semiconductor device capable of estimating fluctuation of the 
circuit characteristics due to fluctuation of the manufacturing process quickly and with high 
precision, so as to suppress fluctuation of the circuit characteristics, which depend on the 
specification of the circuit characteristic. 

Although the embodiment in which the present invention is applied to the MOS 
transistor has been described above, this can be applied to other devices such as Bipolar 
Junction Transistor (BJT) 7 Ferroelectric Random Access Memory (FeRAM). 

Various modifications will become possible for those skilled in the art after receiving 
the teachings of the present disclosure without departing from the scope thereof. 
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